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ABSTRACT 


At NASA, components and subsystems of components in the Shuttle 
and Space Station generally go through a number of redesign stages. 
While data on failures for various design stages are sometimes 
available, the classical procedures for evaluating reliability only 
utilize the failure data on the present design stage of the component 
or subsystem. Often, few or no failures have been recorded on the 
present design stage. 

Last Summer, in the NASA Faculty Fellow Program , Bayesian 
estimators for the reliability of a single component, conditioned on 
the failure data for the present design, were developed. These new 
estimators permit NASA to evaluate the reliability, even when few or 
indeed no failures have been recorded. Point estimates for the later 
evaluation were not possible with the “classical” procedures. 

Since different design stages of a component (or subsystem) 
generally have a good deal in common, the development of new 
statistical procedures for evaluating the reliability, which consider 
the entire failure record for all design stages, has great intuitive 
appeal. 

A typical subsystem consists of a number of different components 
and each component has evolved through a number of redesign 
stages. 

The investigations this Summer considered compound estimation 
procedures and related models. Such models permit the statistical 
consideration of all design stages of each component and thus 
incorporate all the available failure data to obtain estimates for the 
reliability of the present version of the component (or subsystem). 

A number of models were considered to estimate the reliability of a 
component conditioned on its total failure history from two design 
stages. 

It was determined that reliability estimators for the present design 
stage , conditioned on the complete failure history for two design 
stages have lower risk than the corresponding estimators 
conditioned only on the most recent design failure data. 
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Several models were explored and preliminary models involving the 
bivariate Poisson distribution and the Consael Process (a bivariate 
Poisson process) have been developed. Possible shortcomings of the 
models are noted. An example is given to illustrate the procedures. 

These investigations are ongoing with the aim of developing 
estimators that will extend to components (and subsystems) with 
three or more design stages. 
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INTRODUCTION 


Components in the NASA Shuttle and indeed in many other complex 
systems often go through a number of redesign stages. Classical 
reliability estimators rely only on the failure data for the present 
design. Since previous design stages often have a good deal in 

common with the present design, statistical procedures for 
estimating reliability may be improved by also taking into account 
the failure data on the earlier design versions as well. 

Last year in the NASA Summer Faculty Fellow program (SFF), 
Bayesian estimators for the expected value of the reliability of a 
single component, conditioned on its failure history (for the present 
design stage) were developed for the cases of (1) a constant failure 
rate - the exponential model and (2) a variable failure rate - Weibull 
model [1],[2]. 

For the constant failure rate model it was shown that: 

n+1 

E[ R(t)l N(T) = n ] = 1/(1 +t/T) (1) 


where 

R(t) = the reliability or probability that the component will 
successfully function up to time t in the future, 

N(T) = the number of failures up to time T (in the past failure 
history), and 

n = the number of failures recorded up to time T for this 
component. 


These new estimators enable NASA to evaluate reliabilities when 
few or even no failures have been recorded. Evaluation in the latter 
case was not possible with the previous “classical “ estimators. 
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In March 1990, these Bayesian estimators were employed along with 
the classical estimators in a NASA report on the reliability of Orbiter 
APU hydraulic hoses [3]. Both sets of estimators predicted a very 
high reliability for success of the hoses on the next mission (.999..). 

This report will focus on investigations to extend the constant failure 
rate model to utilize the total failure history on a component with 2 
or more design stages. The investigations considered compound 
estimation procedures in order to utilize this failure data history. 
The report incorporates the Bayesian estimators developed last year 
and noted earlier in this report. 

Motivation for the Study 

Since redesign versions of a component would appear to have some 
commonality, the idea of reliability estimators which incorporate the 
total failure history of all design stages of a component seems worth 
considering. When data occurs, for example, on a number of identical 
valves which have been through 2 design stages with the number of 
failures on the j-th valve with the first design denoted by Njj and 

the number on the i-th valve with the new design denoted by 
one can plot Nj versus N 2 The data often suggest some correlation 

between the two failure counts. In general, the number of failures 
Nj and N 2 are not independent. In the discussion to follow, Nj(Tj) 

and N 2 (T 2 > the number of failures up to times Tj and T 2 

recorded for the old and new designs respectively are each assumed 
to have a Poisson distribution with possibly different failure rates j 

and A 2 which are unknown. Since Nj(Tj) and ^^ 2 ) are not 

assumed to be independent, the problem of obtaining Bayesian 
reliability estimators conditioned on the failure data Nj and N 2 

requires some form of a joint or compound probability distribution 
for Nj and N 2 . 

The impact of such estimators that utilize the total failure history 
from 2 or more design stages of a component is indicated by the 
following result. 
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Proposition 


E[R(t)-E(R(t)IN|(T | ), N 2 (T 2 ))]2 S E[R(t)-E(R(t)l N 2 (T 2 ))]2 

i.e., the estimator conditioned on two stages of failure history has a 
lower risk than the corresponding risk for the estimator 
conditioned only on the most recent design failure history. Hence as 
additional failure data on earlier design stages are included, the 
corresponding risk decreases. 


Proof 

E[R(t)-E(R(t)IN,(T,), N 2 (T 2 ))]2 = 


E[R(t)-E(R(t)IN 2 (T 2 ))- (E(R(t)IN 1 (T 1 ),N 2 (T 2 ))-E(R(t)IN 2 (T 2 )))]2 = 

E[R(t)-E(R(t)IN 2 (T 2 ))]2- 2E[(R(t)-E(R(t)IN 2 (T 2 ))) * 

(E(R(t)IN,(T 1 ),N 2 (T 2 ))-E(R(t)IN 2 (T 2 )))] + ( 2 ) 

E[E(R(t)IN 1 (T 1 ),N 2 (T 2 ))-E(R(t)IN 2 (T 2 ))]2 = 

But since 

E[(R(t)-E(R(t)IN 2 (T 2 )))(E(R(t)IN 1 (T 1 ),N 2 (T 2 ))-E(R(t)IN 2 (T 2 )))] = 
E{E[(R(t)-E(R(t)IN 2 (T 2 )))(E(R(t)IN 1 (T 1 ),N 2 (T 2 )) - 
E(R(t)IN 2 (T 2 )))IN j (Tj ),N 2 (T 2 )] } = 

E[E(R(t)IN! (T, ),N 2 (T 2 ))-E(R(t)IN 2 (T 2 ))]2 

then expression (2) becomes 

E[R(t)-E(R(t)INj(T 1 ), N 2 (T 2 ))]2 = 

E[R(t)-E(R(t)IN 2 (T 2 ))]2 - E[E(R(t)IN 1 (T 1 ),N 2 (T 2 ))-E(R(t)IN 2 (T 2 »]2 


Since the last term is non-negative the result is established. Note that 
this result implies that, if the last term is positive, a strict inequality 
holds. 
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DISCUSSION 


The present investigations have considered two preliminary models 
for the distribution of Nj and N2 which lead to reliability estimators. 


Model A : Bivariate Poisson Model 

The first of these employs the bivariate Poisson distribution.. The 
joint probability function for a 2 dimensional Poisson process 
(N 1 (t),N 2 (t)) with cov( Nj(t),N 2 (t)) = “ut is given in [4] by: 


P[N 1 (t) = n 1> N 2 (t) = n 2 ] = f(n 1 ,n 2 ,t) (3) 


where 


min(n 1 , n 2 ) 

f(n 1 , n 2 , t) = exp(-A jt -* 2 t + vX) X ((* j - v) n 1 ' *(A 2 " v )" 2 ’ * * 

j-0 

t n 1 + n 2 -2) (v,) J ) / [(n, -j)! (n 2 - j)! j!] (4). 

Note that in the discussion to follow, in place of t in equation (3) 
one could use T = Ti + J2 which is the total elapsed time for 
failures for both designs. Utilizing expression (4) and Baye’s 
formula, where X 2 (t) = 1 denotes that the component in its second 
design stage is still operating up to time t, one can show that 

E[R(t)IN 1 (T)=n 1 N 2 (T)=n 2 ] = P[X2(t)=l IN^T^m , N 2 (T))=n 2 ] = 


P[N 1 (T)=n 1 N 2 (T)=n 2 IX2(t)=l] P[X2(t)=l]/P[ N^T) = n v N 2 (T) = n 2 ] 

(5) 
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Assuming a joint Bayesian prior distribution for A|,?^ 2 , this expression becomes: 


f S P[Nj(T)=nj N 2 (T)=n 2 IX2(t)=l,A lt A 2 ] P[X2(0=1 1* 2 A 2 ] g(A i * 2 ) d * i, d *2 


J/P[N 1 (T)=n 1 N 2 (T)=n 2 IAi ^ 2 ] g(^i^2^ d?v l , d?v 2 


( 6 ) 


If one assumes the joint prior distribution g(A 1 ^ 2 )to be uniform on finite rectangles with 
independent, i.e. 


(1/0 l) (l/© 2 ) 


for 0 < A 1 < 0 1 and 
0 < ^ 2 <02 


g(A 1> A 2 )= 


0 


elsewhere 


then, 



one can show that the uniform priors, with the limits extended to the 
entire first quadrant, give 


E[R(t)IN 1 (T)=n 1 N 2 (T)=n 2 ] = 


min(n 1 , n 2 ) 

e" vt ^T [1/(1+ t/T) n 2’j +1 ] (vt)Vj! 

| = 0 


min(n^ , n 2 ) 

I (vt) J /j! 

j = 0 


(7) 


Note that in the case when n^=0 failures one has 


E[R(t)IN 1 (T)=0N 2 (T)=n 2 ] = e"' ot ( 1/(1+- t/T) n 2 +1 ) (8) 


where T is the total elapsed time from the start of the first design. 
Notice the similarity of expression (8) to formula (1) developed 
previously for the reliability estimator given just the failure data for 
the second design. 

Also observe that if Nj and N 2 are completely independent then 


E[R(t)IN 1 (T)=n 1 N 2 (T)=n 2 ] = E[R(t)IN 2 (T)=n 2 ] = 1/(1+ t/T) n 2 +1 - 


Thus this model would, in some sense, appear to generalize the 
earlier model(l) which considered only one design stage. 
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Example 


As an application of this reliability estimator consider the follpwing 
example. 

Arbous and Kerrich [5] recorded the number of accidents of 122 
individuals (shunters) in two consecutive time periods. For each 
individual, the number of accidents in the first 6 year period was 
recorded and then after new insurance and safety procedures were 
implemented, the number of accidents for the same individual was 
recorded in the next 5 years. The authors estimated vT ~ .257 for 
the T = 11 years. 

One can use expression (7) to evaluate the “reliability” for the next 
year of an individual, selected at random, given his/her total past 
accident (failure) history. Note that in this case the reliability is just 
the probability that a randomly selected individual with the given 
accident history will not have an accident in the next t years. 

For a randomly selected individual with accident history given by 

m=3 , n2 = 1 , t =1 , -oT ~ .257 , Ti = 6, T2 = 5, T = Ti + T2=l 1 

one finds by using expression (7) that 

E[R(t)IN 1 (ll)=3 N 2 (ll)=l] = .6067624 

Thus the reliability estimator obtained from this model suggests that 
for a randomly selected individual who experienced 3 accidents in 
the first 6 years and 1 in the next 5 years, the probability of no 
accidents in the next year is approximately 61 percent. Such 
information may be used in the setting of insurance premiums for 
the next year for various classes of individuals based on their past 
accident histories. 

Note that in terms of NASA component reliability: 

1. The individuals correspond to different copies of a single 
component with two design stages. 
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2 . Each copy is located in a somewhat similar environment^. g. 
the copies may be located on each of the three space shuttles). 

3 . Data on failures of each copy were recorded in the first design 
stage (i.e. the first 6 years of accident data). 

4 . A second improved design replaced the first and failures were 
recorded. 


Then E[R(t)INj(T^)=ni > N2(T2)=n2 ] gives the probability that the 

second design stage with a given failure history in this location will 
not fail in the next t time units. Note that the NASA description 
assumes that when the component fails it is replaced or repaired so 
that it is equivalent to the original system before the failure. 


It should be pointed out that in Model A some assumptions about the 
relationships between ^ j , X 2 » an< * v were made. In particular, 

v<A j.^2 suggests that v is related to the priors. With this in mind, a 
second preliminary model has been developed. 


Model B : Consael Model 

As was noted earlier, Nj and N2 are probably not independent in 

general. The Consael process [ 6 ] defines a bivariate compound 
Poisson process by 


P[N 1 (T)=n 1 N 2 (T)=n 2 ] = 

JJV*1 T [(A,T) n l]/ni! e-*2T [( * 2T) n2 ]/n 2! * gpi, A^dAjdAj (9) 
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I 


In the Consael process for fixed values of A ^ and A 2, Nj and N2 
correspond to independent Poisson processes while Aj and A 2, have 
a joint density function g(A j A 2). 

If one assumes the triangular density function ( A 2! Aj since one 
assumes that the newer design is an improvement) one has: 


g(A !,A 2 ) = 


2/e j 2 

0 


for 0 < A 1 < 0 j and 
0 < A 2 < A 1 
elsewhere 


( 10 ) 


Utilizing expressions (9) and (10) ,and Baye’s formula in expression 
(6) (and taking limits on the prior distribution) one can again obtain 
an estimate for the reliability: 


E[R(t)IN 1 (T)=n 1 ,N 2 (T)=n 2 ] 


At the present time investigations are continuing with this estimator 
which has a rather complex , highly combinatorial closed form. It is 
anticipated that further consideration of this model will indicate 
approaches to the development of the general model which can 
incorporate failure data from any number of redesign stages . 


SUMMARY AND CONCLUSIONS 


Various NASA failure data suggests support for the development of 
compound/mixture models to estimate the reliability of components 
that have failure data recorded on more than one design stage. 
Bayesian estimators that can utilize all “relevant” failure data, even 
from an earlier design of the component, were investigated. 
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A Bayesian estimator based on the bivariate Poisson distribution was 
developed and an example illustrating the technique applied to 
insurance data was given. The similarity of this example to NASA 
reliability problems was also noted. 

An additional estimator based on the Consael process was developed. 
Investigations into this model are continuing. It is anticipated that 
these investigations will lead to a general model to utilize all relevant 
failure data on a component (or subsystem of components) that has 
experienced more that two design stages. 
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